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ABSTRACT 

We present observations of Can, Znll and Cm absorption lines in 16 damped Lyman 
alpha systems (DLA) and six subDLAs at redshifts 0.6 < z a bs < 1.3, obtained for the 
dual purposes of: (i) clarifying the relationship between DLAs and absorption systems 
selected from their strong Can lines, and (ii) increasing the still limited sample of 
Zn and Cr abundance determinations in this redshift range. We find only partial 
overlap between current samples of intermediate redshift DLAs (which are drawn 
from magnitude limited surveys) and strong Can absorbers: approximately 25 per 
cent of known DLAs at these redshifts have an associated Ca n A3935 line with a rest- 
frame equivalent width greater than 0.35 A, the threshold of the Sloan Digital Sky 
Survey sample assembled by Wild and her collaborators. The lack of the strongest Ca II 
systems (with equivalent widths greater than 0.5 A) is consistent with these authors' 
conclusion that such absorbers are often missed in current DLA surveys because they 
redden and dim the light of the background QSOs. 

We rule out the suggestion that strong Can absorption is associated exclusively 
with the highest column density DLAs. Furthermore, we find no correlation between 
the strength of the Can lines and either the metallicity or degree of depletion of 
refractory elements, although the strongest Can absorber in our sample is also the 
most metal-rich DLA yet discovered, with [Zn/H] ~ solar. We conclude that a complex 
mix of parameters must determine the strengths of the Can lines, including the density 
of particles and ultraviolet photons in the interstellar media of the galaxies hosting 
the DLAs. We find tentative evidence (given the small size of our sample) that strong 
Ca II systems may preferentially sample regions of high gas density, perhaps akin to the 
DLAs exhibiting molecular hydrogen absorption at redshifts z > 2. If this connection is 
confirmed, strong Can absorbers would trace possibly metal-rich, RVbearing columns 
of cool, dense gas at distances up to tens of kpc from normal galaxies. 

Key words: Quasars: absorption lines - galaxies: abundances - galaxies: ISM - 
galaxies: haloes - intergalactic medium 



1 INTRODUCTION 

Quasar absorption lines are a powerful tool for the study 
of gaseous galactic structures over cosmic time, providing 
clues to the nature and evolution of galaxies and galactic 
halos. The large columns of neutral hydrogen gas traced by 
damped Lyman alpha (DLA) and subDLA systems (con- 



* Based on observations made with the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
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ventionally defined to have Hi column densities Af(Hl) ^ 
10 20 - 3 atoms cm" 2 and 10 19 < AT(Hl)< 10 20 ' 3 atoms cm" 2 
respectively) are of particular interest (Rao 2005; Wolfe, 
Gawiser, & Prochaska 2005), as they account for most of 
the neutral gas in the Universe and trace galaxies which are 
often difficult to detect directly. Thus, appropriate to their 
importance, much effort has been expended to understand 
the nature and incidence of the gas traced by DLAs and 
subDLAs. 

While DLAs do contain metals (Wolfe et al. 1986; 
Meyer, Welty, & York 1989; Pettini, Boksenberg, & Hun- 
stead 1990) and ionised gas (Fox et al. 2007a), they are 



2 D. Nestor et al. 



characterized by low metallicities (Kulkarni et al. 2005; Ak- 
erman et al. 2005; Prochaska et al. 2007) and large neutral 
fractions (Viegas 1995; Vladilo et al. 2001). Little, if any, 
evolution is detected in the co-moving mass density of Hi, 
Hhi, from z ~ 4 to z w 0.5 (Rao, Turnshek, & Nestor, 2006; 
Lah et al. 2007; although see Prochaska, Herbert-Fort, & 
Wolfe 2005) and of metals (Kulkarni et al. 2005). This is 
perhaps surprising, considering the putative role of DLAs 
as the reservoir of fuel for star formation. In fact, recent 
work (Wolfe & Chen 2006; Wild, Hewett, & Pettini 2007) 
has shown that in situ star formation in DLAs is signifi- 
cantly less than expected on the basis of the present-day 
Kennicutt- Schmidt law (Kennicutt 1998), considering their 
large Hi surface densities. One possible explanation for the 
apparent lack of evolution in Qui is that a significant frac- 
tion of DLAs with high star formation rates or metallicities 
are missed from traditional magnitude limited optical QSO 
surveys due to chromatic extinction from large columns of 
dust. This possibility, however, is not supported by surveys 
of DLAs in radio-selected QSOs (Ellison et al. 2001; Jorgen- 
son et al. 2006). 

SubDLAs, on the other hand, contribute a smaller frac- 
tion of the total Hi in the Universe. However, they possess 
larger fractions of ionised gas than DLAs, and thus the total 
amount of hydrogen (H I plus H n) in subDLAs may be com- 
parable to that of DLAs (Fox et al. 2007b). Furthermore, 
subDLAs appear to have higher metallicities than DLAs 
(Kulkarni et al. 2007) and thus may be a large repository 
of the metals produced at high redshifts (see Pettini 2006). 
While DLAs and subDLAs both contain large amounts of 
hydrogen, they are likely to arise in different regions within 
galaxies and to trace different environments and properties. 

Strong metal lines of Mg II AA2796, 2804 and 
Fell AA2587, 2600 are ubiquitous in (sub)DLA systems. 
Their large oscillator strengths and near-ultraviolet (UV) 
rest wavelengths make them ideal for tracing large columns 
of neutral gas at redshifts z < 1.7, where the Lya line 
falls below the atmospheric cut-off at 3200 A and therefore 
requires observations from space. In this context, systems 
selected via strong Ca II AA3935, 3970 absorption have re- 
cently garnered attention (Wild & Hewett 2005). Extensive 
literature exists on the Call doublet in the interstellar 
media of the Milky Way and nearby galaxies, since these 
lines can be studied from the ground even at redshift z = 
(e.g. Marshall & Hobbs 1972; Welty, Morton, & Hobbs 1996 
and references therein). This body of work has shown that 
in local interstellar environments Ca is a strongly refractory 
element, exhibiting some of the highest depletion factors 
(Savage & Sembach 1996). Furthermore, since the ionisation 
potential of Ca II is lower than that of hydrogen, Ca II is not 
the dominant ionisation stage of Ca in H I regions, unlike 
Mgn and Fell. Rather, the balance between Call and Cam 
depends on the details of the physical conditions in the gas, 
primarily on temperature and on the densities of particles 
and far-UV photons. Therefore, while all DLAs exhibit 
strong Mgn and Fell lines, the same is not true of Call. 

Wild, Hewett, & Pettini (2006; WHP06) recently con- 
ducted a survey for strong Call systems, with A3935 rest 
equivalent width Wq 3935 ^ 0.35 A, in QSO spectra from the 
Sloan Digital Sky Survey (SDSS). Based on the strengths 
of Znn and Cm lines, WHP06 concluded that, on average, 
strong Can absorbers have iV(Hi) values above the DLA 



limit. They also showed that these absorbers are, on aver- 
age, dustier than typical DLAs, introducing a reddening of 
(E(B ~ V)) > 0.1 magnitudes in the spectra of the back- 
ground QSOs, compared to (E(B - V)) < 0.02 for the gen- 
eral population of SDSS DLAs (WHP06; Vladilo, Prochaska, 
& Wolfe 2008). This level of obscuration of the background 
QSOs is significant, as it causes flux-limited and/or colour- 
selected surveys to underestimate the numbers of strong 
Call systems. If Call-strong DLAs are being missed, and 
if they are preferentially metal-rich relative to DLAs with 
weaker Call absorption, then this class of absorption sys- 
tem may contain a previously overlooked repository of both 
neutral gas and metals. In any case, if Call systems select a 
population of absorbers with unique properties, they could 
hold the key to a better understanding of the detailed phys- 
ical properties of DLAs and subDLAs. 

Since the survey by WHP06, efforts have been devoted 
to better understanding the link between strong Call sys- 
tems and galaxies. Using Jf-band imaging, Hewett & Wild 
(2007) reported an excess of luminous galaxies close to the 
sightlines to 30 QSOs with intervening Call systems at red- 
shifts 0.7 < z < 1.2. The galaxies appear to exhibit a 
luminosity-dependent cross section for Call absorption and 
have a mean impact parameter of ss 25 kpc. Together with 
the incidence of absorbers from WHP06, this result im- 
plies a halo volume filling factor of ~ 10 per cent out to 
~ 35 kpc from the absorbing galaxies. At lower redshifts, 
Zych et al. (2007) were able to associate luminous, metal- 
rich, star-forming spiral galaxies with four out of five strong 
Can absorbers at z < 0.5. 

Although these studies represent progress in uncover- 
ing the nature of Ca II absorbers, ignorance of iV(H i) in such 
systems limits the direct insight they provide into DLAs and 
subDLAs in general. While it seems very likely that Ca II ab- 
sorbers select large columns of neutral gas, direct measure- 
ments of iV(Hl) are currently available for very few Call 
systems. The reason is simple: by the time the SDSS data 
highlighted the potential importance of absorption systems 
selected via Call, no space-borne UV spectrograph was in 
operation to measure the associated column densities of Hi. 
This unfortunate state of affairs is about to change with the 
forthcoming installation of the Cosmic Origins Spectrograph 
(COS) on the refurbished Hubble Space Telescope (HST). In 
the meantime, the most immediate approach to this problem 
is to perform the complementary experiment of measuring 
the strength of Ca II absorption in known DLAs at redshifts 
z £ 1.3, where the Call doublet lines are still accessible with 
optical spectrographs. 

The samples of confirmed DLAs at these intermediate 
redshifts have increased considerably in recent years thanks 
to dedicated HST surveys (see Rao, Turnshek & Nestor 
2006; RTN06). In this paper, we present new observations 
of QSOs from the RTN06 compilation designed to: (i) mea- 
sure the equivalent widths of the Call AA3935, 3970 doublet 
lines in order to assess the overlap of known DLAs (and a few 
subDLAs) with the strong Call absorber sample of WHP06; 
and (ii) simultaneously determine the metallicities and dust 
content of intermediate redshift DLAs via observations of 
the associated Zn II and Cr II absorption lines (Pettini et al. 
1990), thereby adding to the still somewhat limited statistics 
of these measures at z < 1.5 (Akerman et al. 2005; Kulkarni 
et al. 2007). In Section [2] we describe our observations and 
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equivalent width measurements, while Section [3] deals with 
the derivation of column densities and element abundances. 
The main results on the Ca n properties of known DLAs at 
intermediate redshifts are presented in Section [4] We sum- 
marise our principal findings and conclusions in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Sample Selection 

We selected QSOs with known DLAs and subDLAs at 
0.6 < z a b s < 1.3 from the compilation by RTN06. Targets 
were chosen on the basis of absorption redshift (so as to allow 
coverage of the absorption features of interest), and accessi- 
bility during the scheduled observing runs. For each target, 
we first examined the SDSS spectrum at the wavelengths of 
the expected Call absorption. However, those spectra have 
insufficient signal-to- noise ratios (S/N) and resolution to ei- 
ther unambiguously detect or put meaningful limits on Ca II 
absorption strength. 

It is important to remember here that the sample as- 
sembled by RTN06 was constructed from a targeted Lya 
survey of known Mgn absorbers. Therefore, as with any 
absorption line study, the manner in which systems are se- 
lected needs to be considered when interpreting the results 
of our survey. For example, stronger Mg n systems are expo- 
nentially less common than weaker systems (Nestor, Turn- 
shek, & Rao 2005), but more likely to be DLAs (RTN06). 
Combining these two trends, one can demonstrate that the 
unbiased median Wo 2796 for DLAs is ~ 1.35 A, compared to 
1.9 A in our sample. Evidently, our sample is somewhat bi- 
ased against DLAs with weaker Mg n absorption, compared 
to a blind survey. 

Celestial coordinates and magnitudes for each quasar in 
our sample are given in Appendix A. 

2.2 Observations 

Observations were carried out with the ISIS double-beam 
spectrograph on the 4.2 m William Herschel telescope 
(WHT) at the Roque de Los Muchachos Observatory, La 
Palma, in two observing runs of three nights each, in 2006 
April and 2006 November respectively. Atmospheric condi- 
tions varied from partly cloudy to clear, with seeing in the 
range 0.6-2.0 arcsec. In total we secured spectra of 18 QSOs 
covering absorption lines of 16 DLAs and six subDLAs, as 
detailed in Table [T] 

Each observation consisted of a series of simultaneous 
exposures in the blue and red arms of ISIS; individual expo- 
sures typically lasted 1800 s, and the total exposure times in 
each arm ranged from 1 to 5.5 hours per QSO (2-3 hours in 
most cases) . The blue arm was configured so as to cover the 
region encompassing the Zn n and Cr n multiplets, while the 
red arm was centred on the redshifted Call doublet. With 
1200 grooves mm -1 gratings and a 0.8 arcsec entrance slit, 
ISIS delivered a spectral resolution of 0.6-0.7 A full width 
at half maximum (FWHM), sampled with three pixels of 
the EEV (blue arm) and Marconi (red arm) CCDs. Bright 
B-type stars were observed each night to provide templates 
for the removal of telluric absorption lines. The emission 
line spectra of hollow-cathode Cu-Ar and Cu-Ne lamps were 



recorded at hourly intervals to be used for subsequent wave- 
length calibration. Spectrophotometric standards from the 
compilation by Massey et al. (1988) allowed for flux calibra- 
tion of the QSO spectra. 

2.3 Data Reduction and Measurements 

The reduction of the two-dimensional images into co-added 
one-dimensional spectra followed standard procedures. We 
used irai{3 scripts to carry out the initial stages of bias sub- 
traction, fiat-fielding, wavelength calibration, background 
subtraction, and spectral extraction. The individual spec- 
tra of each QSO were co-added with weights proportional 
to the S/N and then binned onto a linear, vacuum helio- 
centric wavelength scale with approximately three bins per 
resolution element. The final steps in the data reduction in- 
volved correction for telluric absorption (if it occurred close 
to absorption lines of interest), and normalisation of each 
spectrum to the local QSO continuum. 

In Figures [T] and [5] we have reproduced normalised 
portions of the spectra encompassing the Znn, Cm and 
Can absorption lines in, respectively, the 16 DLAs and six 
subDLAs, ordered in increasing Right Ascension. As can be 
appreciated from inspection of these two figures, the S/N 
achieved differs significantly from QSO to QSO, depending 
on the variable sky transparency, seeing, QSO magnitude 
and redshifted wavelengths of the features of interest. In ev- 
ery case, however, the combination of S/N and spectral res- 
olution is sufficient to detect a Ca II A3935 absorption line 
with a rest equivalent width Wq 3935 = 0.35 A (the threshold 
of the SDSS survey by WHP06) at the ^ 5<j significance 
level. 

Values of the Wq 3935 are listed in column (5) of Tabled] 
3a upper limits are given for cases where no Ca II absorption 
could be identified. We detected the stronger member of the 
Can doublet, A3935, in 12 out of the 16 DLAs targeted; 
however, only four of the DLAs meet the selection criterion 
Wq 3935 > 0.35 A of the survey by WHP06. We return to 
this point in Section [4] below. For reference, we also give 
in column (4) of Table [T] the rest-frame equivalent width of 
the stronger member of the Mgn doublet, Wo 2796 , from the 
work by RTN06. 

While this study represents the first systematic study 
of Call absorption in (sub) DLAs, UV absorption lines in 
many of the systems in our sample have previously been 
studied by other authors. However, for five of the 16 DLAs 
and one of the six subDLAs, no previous measurements of 
the abundances of Zn and Cr have been reported. For four 
other systems, improved data for Zn and Cr are reported, 
two of which are significant and two which are marginal im- 
provements over published results. Cases where the previous 
measurements are of superior precision to ours are indicated 
by notes in columns (7) and (8) of Table Q] 

As well as considering the individual spectra, we also 
constructed stacked spectra by adding together the nor- 
malised spectra of each DLA in the Znn-Crn and Can 

1 I RAF is distributed by the National Optical Astronomy Observa- 
tories, which are operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 
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Table 1. DLAs OBSERVED, EQUIVALENT WIDTHS, COLUMN DENSITIES, AND ABUNDANCES 
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t from Peroux et al. (2008). Measurements having greater precision are reported by: (a) Peroux et al. (2006a); (b) Meiring et al. 
(2007); (c) Meiring et al. (2006); (d) Meiring et al. (2008); (e) Khare et al. (2004). 



Table 2. Measurements from the literature 



QSO 
(1) 


(2) 


2 abs 

(3) 


W A2796 

(A) 

(4) 


W A3935 

(A) 

(5) 


logiV(Hi) 
(cm -2 ) 
(6) 


logAf(Zn II) 
(cm -2 ) 
(7) 


logiV(Cr II) 
(cm -2 ) 
(8) 


[Zn/H] 

0) 


[Cr/Zn] 
(10) 


Ref. b 
(11) 


0515-4414 


1.713 


1.1508 


2.34 


0.36 ± 0.07 


19.88 + 0.05 


12.18 + 0.02 


12.48 ± 0.04 


-0.33 + 0.05 


-0.72 ± 0.04 


1,2 


0738+313 


0.630 


0.0912 




0.19 + 0.01 


21.18+°,;°, 


< 12.662 


13.28 + 0.22 


< -1.15 


> -0.40 


3 


0738+313 


0.630 


0.2210 


0.61 


0.06 + 0.01 


20.901°;°? 


< 12.825 


13.11 + 0.24 


< -0.71 


> -0.74 


3 


1436-0051 


1.275 


0.7377 


1.11 


0.40 ± 0.02 


20.08l°;i° 


12.67 + 0.05 


< 12.71 


-0.04 + 0.12 


< -0.98 


4 


2328+0022 


1.309 


0.652 


1.90 


0.24 + 0.02 


on qo+0.06 


12.43 + 0.15 


13.35 + 0.19 


-0.52 + 0.22 


-0.10 + 0.24 


5 


2331+0038 


1.486 


1.1414 


2.53 


0.26 + 0.02 


20.00t°;° 4 


12.22 + 0.09 


< 12.37 


-0.41 + 0.12 


< -0.87 


6 


2335+1501 


0.790 


0.6798 


0.97 


0.19 + 0.03 


I9.70l°; 3 3 ° 


12.37 + 0.04 


12.89 + 0.10 


+0.04 + 0.30 


-0.50 + 0.11 


7 



b References— 1: de la Varga et al. (2000); 2: Quast et al. (2008); 3: Meiring et al. (2006); 4: Meiring et al. (2008); 5: Peroux et al. (2006a); 
6: Meiring et al. (2007); 7: Peroux et al. (2008). 



regions (in the DLA rest frame), following the same pro- 
cedure as WHP06. The resulting average spectrum, which 
includes DLAs in which the Call lines are below our de- 
tection limit, is useful for comparative studies, as discussed 
below. A second average spectrum was assembled by adding 
together the data for each subDLA. The borderline DLA 
towards 1107+0003 (log N (Hi) = 20.26 ± 0.1) was included 
in the DLA stack; the inclusion/exclusion of this system 
in the DLA/subDLA stack has little effect on the results. 
The Znll-Crn region of the DLA towards 0449-168 was 
excluded from all stacks, as it falls within the Lya forest 
(see Figure [TJ. 



3 COLUMN DENSITIES AND ELEMENT 
ABUNDANCES 

Column densities of Zn II and Cr II were measured by fitting 
the observed absorption lines with theoretical profiles gen- 
erated with the software package VPFIT0 as described, for 
example, by Rix et al. (2007). Briefly, vpfit uses x 2 min- 
imisation to deduce the values of redshift z, column density 
N (cm -2 ), and Doppler parameter b (km s" 1 ) that best re- 
produce the observed absorption line profiles, vpfit takes 



2 vpfit is available from |http: //www. ast.cam.ac.uk/ 
rfc/vpfit.html 
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Figure 1. The Znll-Crn (left) and Can (right) regions of the spectra of the 16 DLAs observed. The spectra have been normalised 
to the underlying QSO continuum and reduced to the rest frame of the DLA. In each panel, the horizontal dash-dot line indicates the 
continuum level, while the dotted line near the bottom of each plot shows the lc error spectrum. Vertical long-dash lines are drawn at 
the wavelengths of the Zn IIAA2026, 2063 doublet, Cr II AA2056, 2062, 2066 triplet, and Ca II AA3935, 3970 doublet. 



into account the instrumental broadening function in its 
X 2 minimisation and error evaluation. Laboratory (vacuum) 
wavelengths and /-values of the transitions analysed were 
taken from the compilation by Morton (2003) with subse- 
quent updates by Jenkins & Tripp (2006). In most cases, the 
resolution of our spectra (FWHM w 30-60 km s" 1 ) is insuf- 
ficient to resolve the intrinsic velocity structure of the gas. 
However, as we are dealing with weak absorption lines, we 
do not expect line saturation to be a problem (see the dis- 
cussion by Pettini et al. 1990); the measured doublet (Znn) 
and triplet (Cr n) ratios confirm that this is indeed the case 
for the systems analysed here. 

Values of TV(Znn) and TV(Crll) are listed in Table [H 
together with the errors estimated by vpfit on the basis 
of the S/N ratios of the spectra. Since these species are 
the dominant ionisation stages of Zn and Cr in H I regions 



(e.g. Morton et al. 1973; Viegas 1995), the corresponding 
element abundances can be inferred directly from the ra- 
tios iV(Znn)/JV(Hi) and JV(Grii)/JV(Hi). Values of AT (Hi) 
for the systems observed are reproduced in column (6) of 
Table [1] from the compilation by RTN06. Reference to the 
solar values of Zn/H and Cr/Zn (log (Zn/H) Q = —7.37 and 
log (Cr/Zn) Q = +1.02; Lodders 2003), then gives the quan- 
tities [Zn/H] and [Cr/Zn] in the usual definition [X/Y] = 
log (X/Y) DLA — log (X/Y) Q . The former is often adopted as 
a measure of the metallicity, while the latter provides an 
indication of the degree of dust depletion of refractory el- 
ements, for the reasons discussed by Pettini et al. (1990). 
Values of [Zn/H] and [Cr/Zn] for the systems observed are 
given in the last two columns of Table [1] the errors quoted 
were obtained by combining in quadrature the errors on the 
column densities of the respective ions. 
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As noted above, measurements of Zn and Cr abun- 
dances have been previously reported in the literature for 
many of the absorbers in our sample. Out of the ten DLAs 
in common with previous work for which we have esti- 
mates of [Zn/H] and [Cr/Zn], eight are in very good agree- 
ment with published values. The exceptions are the DLAs 
in the spectra of the QSOs 1225+0035 and 1733+5533 (at 
Zabs = 0.7728 and 0.9981 respectively) where the values of 
[Zn/H] implied from our detections of the Znll doublet lines 
(see Figure [TJ are actually higher than the upper limits es- 
timated by Meiring et al. (2006). However, these authors' 
data are of rather poor S/N at the relevant wavelengths. 



3.1 Metallicity Evolution of DLAs 

The most recent compilation of [Zn/H] measurements in 
DLAs by Kulkarni et al. (2007) includes 18 systems in the 
redshift range 0.6 < z < 1.3 (corresponding to a time inter- 



val amounting to nearly one quarter of the current age of 
the universe). Thus, the five new (and two improved) mea- 
surements reported here constitute a significant addition to 
the DLA statistics at these epochs. In Figure [3] we plot the 
full sample of [Zn/H] measurements as a function of red- 
shift. The new/improved data reported here, shown with 
solid squares, confirm the wide range of metallicities exhib- 
ited by DLAs at all redshifts: as pointed out by Pettini et 
al. (1994, 1997b), at any one epoch DLA metallicities can 
span two orders of magnitude, from nearly solar to ~ 1/100 
of solar. 

Also shown in Figure [3] are the mean values of 
the column density-weighted metallicity, [(Zn/H)] (where 
(Zn/H) = J2 N (. Znll )/J2 N(Ri)— see Pettini et al. 1994), 
for the redshift range of our sample (0.6 < z a b s < 1-3) and 
for higher redshift DLAs in two bins: 1.3 < z a b s < 2.6 and 
2.6 < z a bs < 4. As discussed by Kulkarni et al. (2007 and 
references therein), there is evidence for mild evolution of 
the mean DLA metallicity with the progress of time, in the 
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Figure 2. As in Figure [T] but for the six subDLAs observed. 



sense that the value [(Zn/H)] = —0.71 ± 0.08 we measure 
at 0.6 < Zabs St 1.3 is a factor of about two to three higher 
than in the other two redshift bins: [(Zn/H)] = -1.06 ±0.06 
at 1.3 < Zabs «S 2.6, and [(Zn/H)] = -1.14 ± 0.11 at 
2.6 < Zabs < 4.0, respectively. 

One of the factors which affect this type of compari- 
son, and yet is generally overlooked, is the fact that DLAs 
with z a b s < 1.7 are normally selected via their (strong) 
MgnA2796 absorption, while the selection at higher red- 
shifts is via the damped Lya line itself. The two criteria 
may not be equivalent; in particular, it now seems well es- 
tablished that DLAs with strong Mgn absorption are on 
the whole more metal-rich than the average (Nestor et al. 
2003; Murphy et al. 2007). To assess the importance of this 
potential bias, we use four distributions: (i) the distribution 
of DLA metallicities as a function of Wq 2796 - at redshifts 
0.6 < z < 1.3, we find that [Zn/H] oc (Wo 2796 ) 1 7±0 1 is a 
good representation of the 12 measured values in our sample 
(Pearson's r = 0.85 ± 0.08, \l = 2.7); (ii) the DLA fraction 
as a function of Wq 2796 from RTN06; (Hi) the distribution 
of values of Wq 2796 in SDSS QSOs from Nestor et al. (2005); 
and (iv) the distribution of values of Wq 2796 in the sample 
used to estimate [(Zn/H)], to arrive at the conclusion that 
the value of [(Zn/H)] = -0.7 ± 0.08 at 0.6 < z a bs < 1.3 is 
probably biased too high by ~ 0.1 dex. Although this is not 
a large factor, it nevertheless dilutes further the evidence for 
redshift evolution of the mean metallicity of DLAs. 



3.2 Depletion of Chromium 

In Figure [4] we have plotted [Cr/Zn] as a function of DLA 
metallicity (as measured by the [Zn/H] ratio), again combin- 
ing the systems in the compilation by Kulkarni et al. (2007) 
with the new, or improved, measurements from the present 
work. Both Zn and Cr track the abundance of Fe in Galactic 
stars with [Fe/H] ^ — 2 (i.e. in the same metallicity regime 
as the vast majority of DLAs — see Cayrel et al. 2004 and 
Nissen et al. 2007) . On the other hand, in the ISM near the 
Sun, Cr is among the most depleted elements with only a 
small fraction in the gas phase while most of it is 'hidden 
away' in dust grains; Zn, on the contrary, exhibits at most 
only mild depletions (Savage & Sembach 1996). Thus, the 
[Cr/Zn] ratio in DLAs is a convenient way to assess the de- 
gree to which refractory elements are depleted in the ISM 
of the galaxies hosting these absorbers (Pettini et al. 1990) 
and, by inference, obtain an indication of their dust content 
(Pettini et al. 1997a). 

In Figure |3] we see the clear correlation between Cr de- 
pletion and metallicity that has already been noted (e.g. Pet- 
tini et al. 1994; Ledoux, Petitjean & Srianand 2003; Aker- 
man et al. 2005; Meiring et al. 2006); the data reported here 
fit the known trend and improve the statistics at the high 
metallicity end. Although the measurements scatter about 
the mean relation between Cr depletion and Zn abundance, 
the data in Figure [4] show that in DLAs with metallicities 
[Zn/H] ^ — 1.5, Cr suffers little depletion relative to Zn; 
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Figure 3. Measurements of [Zn/H] in DLAs at < z a b s < 4 
from the compilation by Kulkarni et al. (2007) and from the 
present work. Except for the cases in which there exists a su- 
perior measurement in the literature, the DLAs from the present 
study are shown as red squares. Upper-limits are indicated with 
downward pointing arrows. Also shown are the mean values of 
column density-weighted metallicity [(Zn/H)] (see text for defi- 
nition) in three redshift intervals, together with the ±1<t ranges 
calculated via bootstrap techniques. As discussed in the text, the 
value of [(Zn/H)] at 0.6 < 2 a b s ^ 1.3 may be biased too high 
by about 0.1 dex, because these DLAs were primarily selected via 
their strong Mgn absorption. The vertical bar in the lower left- 
hand corner gives an indication of the typical error in [Zn/H]. 





1 1 1 1 1 1 1 

x t 


X 


4 




*;* 4 x >< >< : 




X X x< x 




x x 

X>< x X >< 

x x 




m x 
x x X x "x X ■ 




x ■ 

* ■ 




X x 

X x ■ x x 

1 X 




■ X 



-2 -1.5 -1 -0.5 



[Zn/H] 

Figure 4. Measurements of [Cr/Zn] in DLAs from the compila- 
tion by Kulkarni et al. (2007) and from the present work. DLAs 
from the present work are shown by red squares unless a supe- 
rior measurement exists in the literature. The error bars in the 
lower left-hand corner show the typical uncertainties in the ratios 
plotted. 

in the gas phase (although this would be in contradiction to 
the finding that strong Call absorbers induce higher-than- 
average reddening in the background quasar). We now con- 
sider these possibilities in turn, in light of the new data 
presented here. 



on the other hand, in the DLAs with near-solar metallic- 
ities, ~ 75 per cent of the Cr is presumably incorporated 
into dust ([Cr/Zn] ~ —0.6). This is still a relatively mod- 
est degree of depletion compared to that exhibited by Cr 
in diffuse clouds in the Milky Way disk, and more akin to 
the values seen along sight-lines which sample mostly the 
interstellar medium in the halo of our Galaxy (see Figure 
6 of Savage & Sembach 1996). The origin and implications 
of this metallicity-dependent depletion pattern have been 
considered by Vladilo (2002). 



4 CA II ABSORPTION IN DLAS 

A primary aim of the observations presented in this paper is 
to establish the Ca II content of known DLAs at intermedi- 
ate redshifts. Ultimately, we hope to clarify the relationship 
between absorbers with strong Call lines, which can be se- 
lected from ground-based QSO spectra such as those from 
the SDSS without reference to the Lya line, and the more 
general population of DLAs. 

For the reasons discussed by WHP06, we suspect that 
strong Ca II absorbers may be a subset of the DLA popula- 
tion. However, the properties of DLAs that lead to strong 
Ca II absorption are not clear. WHP06 considered three pos- 
sibilities; they proposed that strong Call absorbers select 
DLAs with: (i) the largest iV(Hl) values; (ii) the highest 
metallicities; or (iii) the largest volume densities, 7i(Htot). 
A further possibility is an absence or destruction of dust 
grains, leading to unexpectedly large fractions of Ca atoms 



4.1 Call as a Tracer of Column Density? 

We examine in Figure [5]the relationship between the equiv- 
alent width of the stronger member of the Call doublet, 
Wq 3935 , and the column density of neutral hydrogen iV(Hl). 
We have augmented our sample with a trawl of the litera- 
ture for other measurements of Call in DLAs and subDLAs, 
collected in Table [2] However, it must be borne in mind that 
the six additional systems thus found do not represent a sta- 
tistical sample, because non-detections of Ca II and the cor- 
responding upper limits on Wj 3935 generally go unreported 
in the literature. 

The data collected in Figure [5] do not support the 
hypothesis that strong Call absorption corresponds exclu- 
sively to DLAs at the upper end of the A r (Hl) distribution. 
WHP06 showed that, if this were the case, the incidence of 
strong Ca II absorbers per unit redshift would match that of 
DLAs with 10 21 < 7V(H l)< 10 22 so that, in this scenario, the 
former would account for all of the latter. In contrast, our 
combined sample includes eight DLAs with iV(H i) ^ 10 21 , 
and none of them is a strong Ca II absorber. 

More generally, inspection of Figure [S] shows that there 
is no simple relation between Wq 3935 and iV(Hl) in DLAs. 
Selecting absorption systems by Call equivalent width does, 
however, primarily select DLAs: among the systems we have 
observed, all those with Wq 3935 > 0.15 A are DLAs [within 
the accuracy of the determination of iV(Hl)]. On the other 
hand, there are subDLAs with Wq 3935 > 0.15 A reported 
in the literature. Conversely, selecting by iV(Hl) does not 
ensure that a strong Call line is present: over one third 
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Figure 5. The rest-frame equivalent width of the CaoA3935 
line plotted as a function of the column density of neutral hy- 
drogen. Open circles: DLAs from the present survey; the size of 
the circle is proportional to the equivalent width of Mg II A2796 
(0.7 A < W^ 2796 < 3.1 A— see Tabled}. Upper limits to W^ 3935 
are 3a. Crosses: additional measurements from the literature 
(listed in Table [S). Note, however, that these are likely to be bi- 
ased towards higher-than-average values of Wq 3935 . Filled (blue) 
triangle and filled (green) diamond: values measured from the 
stacked spectra of, respectively, the DLAs and sub-DLAs in Ta- 
ble [l] plotted at the corresponding mean values of TV (Hi) (see 
section 12.31 1 . The vertical dash line shows the conventional dis- 
tinction between DLAs and subDLAs at AT (Hi) = 2 X 10 20 cm -2 . 



three quarters have Wq 3935 < 0.3 A. On the basis of our 
composite spectrum (blue triangle in Figure[5] — see sec- 
tion 2.3), we conclude that Call absorption with an average 
(Wo 3935 ) = 0.26 A is associated with known DLAs. Fur- 
thermore, despite the few strong detections reported in the 
literature, the composite spectrum of the sub-DLAs in our 
sample shows that their associated Ca n absorption is weak, 
with an average (Wq 3935 ) < 0.02 A (green diamond in Fig- 
ure [5]). 

We note that only four out of the 19 DLAs in the com- 
bined sample meet the selection criterion Wq 3935 ^ 0.35 A of 
the SDSS survey by WHP06, and only one out of 19 would 
be classed as a 'strong' Call absorber (Wq 3935 > 0.5 A) by 
WHP06. Evidently, there is little overlap as yet between con- 
firmed DLAs and the SDSS candidates identified via Call. 
This is not surprising, however, when we consider the mag- 
nitude limit of the HST observations which RTN06 used to 
assemble the current sample of intermediate redshift DLAs. 

Specifically, WHP06 estimate that the true (unbiased) 
redshift path density of Can with Wq 3935 > 0.5 A is ~ 20- 
30 per cent of that of DLAs. However, for a flux-limited 
sample in the SDSS i-band they also calculate that ~ 40 
per cent of their strong Call systems are missed due to 
the relatively strong extinction of the background QSOs 
by such absorbers. A detailed description of the calcula- 
tion is provided in Section 6.1 of WHP06. The majority of 
quasars included in the RTN06 compilation were discovered 
either using much bluer passbands (usually B or V), or from 
an SDSS <?-band flux-limited sample. The fractional loss of 



strong Call-systems from the RTN06 compilation may be 
calculated using exactly the same procedure as described in 
WHP06 but using observed-frame wavelengths appropriate 
for the B-, V- or 3-bands instead of the ~ 8000 A SDSS 
i-band. The significantly shorter absorber rest-frame wave- 
lengths probed by the bluer passbands mean that the frac- 
tion of background QSOs falling below a pre-defined flux- 
limit increases to ~ 80 per cent, compared to the ~ 40 
per cent loss for the SDSS i-band. Assuming that A/(Hi) 
and Wq 3935 are unrelated in DLAs, we thus calculate that 
our sample of 16 DLAs was likely to include only one Can- 
absorber with Wq 3935 > 0.5 A, as indeed turned out to be 
the case (the Wq 3935 = 0.82 A absorber at z abs = 0.7156 to- 
wards 1323 — 0021). Estimates of the degree of reddening for 
individual absorbers in our sample are given in Appendix A. 

Returning to Figure [5] it can also be appreciated that 
there is no obvious correlation between Wq 3935 and Mgn 
Wo 2796 - the size of the open circles in the Figure is pro- 
portional to the latter quantity. Mg 11 absorption in DLAs is 
typically at least moderately saturated and, thus, primarily 
a tracer of the velocity structure of the absorbing gas. 



4.2 Call as a Tracer of Metallicity or Dust 
Content? 

We have searched for other clues to the origin of strong 
Can absorbers by testing for correlations of Wq 3935 with 
DLA metallicity (as measured by [Zn/H]) and dust deple- 
tions (via [Cr/Zn], as explained in Section 13.21 above). As 
can be seen from Figure [6] no such correlation with [Zn/H] 
is evident in the present data set. The only DLA in our 
sample with Wo 3935 > 0.5 A is indeed the most metal-rich 
known at any redshif10, with [Zn/H] = +0.12T0.26 (see Fig- 
ures [6] and . However, DLAs with Wq 3935 > 0.15 A span 
essentially the full range of values of [Zn/H]. Due to the 
dearth of strong Can absorbers in our sample, it is difficult 
to test the hypothesis that the strongest Ca 11 absorbers rep- 
resent the DLAs with the highest metallicities. However, the 
highly-enriched nature of the lone strong Can absorber in 
our survey suggests that metallicity may play some role. 

Given the highly refractory nature of Ca, we also con- 
sider whether a lack of dust grains, leading to a large frac- 
tion of Ca remaining in the gas phase, may be an important 
factor in determining Can absorption strength. However, 
Figure [7] shows no obvious relationship between Can ab- 
sorption and depletion of Cr present. The values of [Cr/Zn] 
determined here are similar to those of other DLA sam- 
ples and, in particular, we find no evidence for higher values 
of Wq 3935 in DLAs with low dust depletions (as evidenced 
by near-solar values of [Cr/Zn]). The lone strong Can ab- 
sorber in our sample has a [Cr/Zn] value consistent with 
the average for systems of similar Ca n-strength reported by 
WHP06. 



3 This system is discussed by Peroux et al. (2006b) as a subDLA 
based on the 7V(H 1) measurement of Khare et al. (2004). However, 
we use the subsequent measurement by Rao, Turnshek and Nestor 
(2006) of AT(Hl)= 20.54 ±0.15. 
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4.3 Call as a Tracer of Volume Density? 

The strengths of Mgn, Znll, and Cm lines in DLAs are 
determined primarily by the H I column density, metallicity, 
and dust depletion, as well as the velocity dispersion of the 
gas and the oscillator strengths of the transitions. This is be- 
cause these ions are the dominant species of Mg, Zn, and Cr 
in Hi regions, where photons with sufficient energy to ionise 
these elements beyond their first ions do not penetrate. The 
same is not true of Call. With an ionisation potential of 
11.87 eV, a significant fraction of Ca can be doubly ionised 
even when the bulk of the gas is neutral. In photoionisation 
equilibrium, the balance between doubly and singly ionised 
Ca is given by: 



i(Ca ; 



2+^ 



n(Ca+) 



*(T)n(e) 



(1) 



where n denotes volume densities, V is the Ca + photoionisa- 
tion rate, ct(T) is the recombination coefficient to form Ca + , 
and n(e) is the electron density. Presumably, this depen- 
dence of n(Ca + )/n(CaTOT) on the local physical conditions 
is the extra factor responsible for the lack of any obvious 
correlation of Wq 3935 with Hi column density, metallicity 
and depletion levels. 

We can exploit this feature of Call absorbers to gain 
information on the physical conditions which determine 
the ionisation balance of Ca. Since we know, or have es- 
timates of, AT (Hi), metallicity (from [Zn/H]), and dust de- 
pletion (from [Cr/Zn]), we can use the measured values of 
JV(Call) - determined via Wo 3935 - to estimate the ratio 
n(Ca 2+ )/n(Ca + ). For a range of values of the intensity of 
the far-UV radiation field (which determines the photoion- 
isation rate F) we can then deduce the corresponding val- 
ues of gas density 72(Htot) (on which n(e) depends), and 
thereby assess the hypothesis the gas volume density is the 
important factor regulating the strength of Call absorption. 
This, in turn, may shed light on the nature of the strong Ca II 
absorbers and their relationship, for example, to DLAs with 
associated molecular hydrogen absorption (Srianand et al. 
2005). 



4-3.1 Photoionisation models 

With this aim, we ran series of photoionisation models with 
the software package CLOumQ (Ferland et al. 1998; Ferland 
2000), approximating the DLAs as slabs of constant den- 
sity gas irradiated by the metagalactic ionising background 
(Haardt & Madau 2001) and the cosmic microwave back- 
ground at the appropriate redshift. We further added cosmic 
rays and an interstellar radiation field due to the OB stars 
within the galaxies giving rise to the DLAs; the spectral 
distribution of this last component was assumed to be the 
same as that of the interstellar radiation field in the solar 
vicinity, as given by Black (1987). For each DLA in Tables [T] 
and [2] where Znii and Cm lines were detected, we adopted 
the metallicity implied by the measures of [Zn/H] (thereby 
assuming the dust depletion of Zn to be negligible in DLAs) 
and included dust grains with an abundance scaled with the 
DLA metallicity. We explicitly accounted for the depletion 
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Figure 6. The rest-frame equivalent width of the Call A3935 line 
plotted as a function of DLA metallicity. The symbols have the 
same meaning as in Figure \5\ except that the size of the open 
circles is now proportional to Af(Hl). SubDLAs are shown in 
green; for some of these, the assumption N(Zn u)/N (H i) = Zn/H 
adopted in section [3] may not be correct, possibly leading to an 
overestimate of the metallicity if some of the Znll is associated 
with ionised gas. 



of metals onto grains by assuming solar relative abundances 
as given by Lodders (2003), and a depletion pattern moti- 
vated by the well-studied interstellar cloud in front of the 
star £ Oph (Savage & Sembach 1996). We scaled the level 
of depletion for each element more highly refactory than Zn 
based on our [Cr/Zn] measurement: 



[X/Zn]= A(X) x [Cr/Zn] , 



(2) 



where A(X) is determined from the depletion values given 
by Savage & Sembach. This functional form matches the 
relative degrees of depletion from the gas-phase in the ISM 
towards £ Oph, while reducing to [X/Zn] gas = when 
[Cr/Zn] gas = (i.e. when Cr is undepleted, so are all other 
metals). 

Each model was characterised by a pair of values: gas 
density and intensity of the external stellar radiation field. 
We then used cloudy to predict the value of A^Caii) ap- 
propriate to that model given the measured A^Hl) and all 
of the above parameters. We assigned a \ 2 confidence level 
to each model, by comparing predicted and observed val- 
ues of A r (Call) and taking into account the uncertainties in 
our measurements of column densities. We assumed an un- 
certainty of ±0.3 applied to the scaling factor in eq.© for 
Ca; i.e., A(Ca) — 1.9 ± 0.3. Although ad hoc, it is likely a 
conservative estimate. The relative contribution to the total 
error arising from this uncertainty depends on the level of 
depletion, and ranged between zero and 50 per cent, with an 
average of 25 per cent, in our models. Thus, for each DLA in 
Tables 1 and 2 with measured A r (Zn n) and A^Cr n) we were 
able to produce confidence intervals in gas density (cm -3 ) 
for a given choice of intensity of the interstellar radiation 
field (in units of the local ISM radiation field intensity) . 

Figure [8] shows the la ranges in gas density if the in- 
tensity of the interstellar radiation field is allowed to vary 
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Figure 7. The rest-frame equivalent width of the Can A3935 line 
plotted as a function of Cr depletion. The symbols have the same 
meaning as in Figure [6] 



between 1/10 and three times the reference value appropri- 
ate to our location within the Milky Way. We see that, even 
in the best cases, we only constrain n(HxoT) to within « 1-2 
orders of magnitude, depending on the level of star forma- 
tion activity in the galaxies hosting the DLAs (since the star 
formation rate is presumably the main determinant of the 
radiation field intensity). Previous estimates of the interstel- 
lar radiation field intensity in DLAs have been deduced from 
consideration of the relative populations of the rotationally 
excited levels of the H2 molecule (e.g. Ge & Bechtold 1997). 
Available data indicate a UV field intensity comparable to 
the local value (Srianand et al. 2005), consistent with the 
conclusion from most imaging studies that the galaxies giv- 
ing rise to DLAs are generally fainter than L* (Rao et al. 
2003 and reference therein). 

These measurements refer to DLAs at redshifts z ^ 2 
and it is thus unclear to what extent they also apply to lower 
redshift systems. If we assume, as a working assumption, a 
UV field equal to that observed locally, we obtain the values 
of gas density shown in Figure [9] and listed in Table [3] (where 
the DLAs are now listed in increasing order of Wo 3935 ). The 
range of Wq 3935 in our sample is too small, and in some cases 
w(Htot) is too poorly constrained, to discern an overall cor- 
relation between Wq 3935 and n(HTOT)- However, we note 
that in the DLA with the lowest measured value of Wq 3935 
(at z abs = 0.8591 in 0454+039, with Wq 3935 = 0.063 A) we 
deduce the lowest value of gas density, n(HxoT) — 0.1 cm -3 . 
Such a weak Can line cannot be ascribed to dust depletion 
of Ca, since in this DLA the Cr/Zn ratio is solar. Conversely, 
the only 'strong' Call absorber in our sample, according to 
the definition by WHP06 (at z abs = 0.7156 in 1323-0021, 
with Wq 3935 = 0.82 A), requires relatively high density ac- 
cording to our modelling, with n(HTOT) — 7cm -3 . In this, 
as in other DLAs with large values of Wo 3935 , we may have 
underestimated iV(Call), and therefore n(HTOT) (by ~ 0.1 
- 0.5 dex), if there is hidden saturation in the absorption 
lines. 

We can apply a similar analysis to the 'low- Wq 3935 ' 
«W£ 3935 ) = 0.49 A) and 'high-W^ 3935 ' ((W^ 3935 ) = 1.01 A) 




-2 -1 0.1 2 

Log[gas density (cm )] 

Figure 8. The rest-frame equivalent width of the Call A3935 line 
plotted as a function of the values of n(HxoT) deduced from our 
CLOUDY modelling. The range of values of ii(Htot) shown cor- 
responds to the full range of normalisations of the UV radiation 
field considered in our modelling, from 1/10 to three times the lo- 
cal interstellar value. Upper-limits for Wq 3935 and n(HxoT) are 
shown at one-sigma. 



-1 

Loglgas density (cm 



')] 



Figure 9. As for Figure[8] but now assuming that the interstellar 
radiation fields in the galaxies hosting the DLAs have the same 
intensity as that of the Milky Way near the Sun. 



subsamples of Call absorbers from WHP06, for which those 
authors measured values of A r (Call), iV(Znll), and iV(Crn) 
from composites of SDSS spectra. Since we do not know the 
value of Af(Hl) appropriate to these composites, we have 
performed the calculation for a range of values of the neu- 
tral hydrogen column density (and therefore also a range 
of metallicities) . For typical DLA A jr (H l)-values, we find 
that, for a normalisation of the UV field equal to the lo- 
cal value, the typical strong Call absorber arises in gas with 



density in the range n(HxoT) 



10 u 



10 1 ' 3 atoms cm 
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for the low-W^ 3935 subsample, and ti(Htot) ^ 10 1 1 - 
10 1 ' 6 atoms cm -3 , for the high-Wo 3935 subsample, depend- 
ing on the value of JV(Hi). 

4-3.2 Connection with H2? 

The current evidence thus suggests that strong Call ab- 
sorbers trace gas of relatively high density. Returning to our 
sample, it can be seen from Figure [9] and Table [3] that if 
the UV radiation field intensity is comparable to that in 
the solar neighbourhood, most Can-bearing DLAs, and all 
four with W^ 3935 > 0.35 A, likely arise in gas with densities 
greater than ~ 1 atom cm -3 , while for the average Can 
system of WHP06 we find densities > 10 atom cm -3 . These 
values approach the range derived by Srianand et al. (2005) 
in DLAs at z > 2 also exhibiting H2 absorption, from con- 
sideration of the relative populations of the fine-structure 
levels of the ground state of C I. These authors propose that 
such systems trace the 'Cold Neutral' phase of the inter- 
stellar medium in galaxies; our study raises the possibil- 
ity that absorption systems selected from their strong Call 
absorption may be the counterpart of such DLAs at lower 
redshift. In the last column of Table [3] we list the molec- 
ular fraction, /h 2 = 2N(H 2 )/(2 N(H 2 )+iV(H 1)), predicted 
by our CLOUDY models for a normalization of the local 
UV field of unity. The range of values is broadly consistent 
with those found by Noterdaeme et al. (2008) for DLAs at 
redshifts 2 < z < 4, although we predict a smaller fraction 
of very low /H 2 -values. However, the lowest /H 2 -values in 
the Noterdaeme et al. sample are dominated by the systems 
with the lowest metallicity and depletion, both of which are 
(on average) greater at the redshifts of our sample (Fig- 
ures [3] and |3| . If we consider only systems with metallici- 
ties greater than ~ —1.3 in the high-z Noterdaeme et al. 
sample, the values predicted by CLOUDY are in reason- 
able agreement with observed molecular fractions. For DLAs 
with properties similar to the Iow-Wq 3935 and high-Wo 3935 
subsamples of WHP06, our models predict /h 2 — —2.8 to 
—2.1 and /h 2 — —0.4, respectively. If these interpretations 
are correct, we expect that future space-based observations 
of strong Call absorbers at far-UV wavelengths will show 
not only a damped Lya line, but also detectable levels of H2 
absorption. 

At present, we can test this prediction with only three 
reported measurements of, or limits on, iV(H2) at z a b s < 1.3: 
in the DLAs towards 0454+039 (which is in our sample) 
and 1328+307 (z abs = 0.692; AT (Hi) = 10 21 ' 28 cm" 2 ), and 
in the subDLA towards 0515-4414 (z abs = 1.1508; iV(Hl) = 
10 19 ' 88 cm -2 - see Table 2.) Molecular hydrogen is unde- 
tected in either of the two DLAs (Ge & Bechtold 1999), 
with N(U 2 ) < 10 16 - 38 and /h 2 < 10~ 4 01 for the system to- 
wards 0454+039 and N(H 2 ) < 10 15 ™ and fn 2 < 10" 5 28 for 
the system towards 1328 + 307. The former has the weakest 
detected Call line in our sample. For the latter, we esti- 
mate Wo 3935 < 0.2 A from its SDSS spectrum. Thus, nei- 
ther of these two DLAs are strong Can absorbers. In con- 
trast, the subDLA towards 0515-4414 does exhibit H2 ab- 
sorption (Reimers et al. 2003), with JV(H 2 ) = 10 16 - 94±0M 
and /h 2 = l(T 2 ' 64±0 - 30 . In this system, Wq 3935 = 0.36 A; 
while this does not qualify it as a 'strong' Call absorber 
in the sense defined by WHP06, we note that this value of 
Wo 3935 is unusually high considering the low iV(Hl) (see 



Table 3. Predicted gas densities and molecular fractions' 1 



QSO 


■\yA3935 

(°A) 


log niH^vci^r ) 
(cm -3 ) 




3 J £i '2 


0449- 168 b 


< 


0.03 


< -1.6 


< 


-9.4 


1733+5533 b 


< 


0.02 


< -0.45 


< 


-6.6 


1007+0042 b 


< 


0.03 


< -0.14 


< 


-5.5 


0454+039 


0.06 


±0.01 


1 - ua -0.18 




-8.5 


1727+5302 


0.16 


±0.02 


0.71 ±0.18 




-3.9 


2353-0028 


0.17 


±0.03 






-1.2 


2328+0022 


0.24 


±0.02 


-°-ob12:E 




-6.3 


1010+0003 


0.26 


±0.02 


0.47 ±0.13 




-3.6 


1727+5302 


0.29 


±0.03 


0.22 ±0.20 




-3.1 


1107+0048 


0.34 


±0.01 


0.31 ±0.21 




-3.1 


1225+0035 


0.36 


±0.03 


n QQ+0.22 
u.oa_ Q 2g 




-2.6 


0256+0110 c 


0.39 


±0.04 






-3.3 


1323-0021 


0.82 


±0.02 


0-82«; 3 i 




-0.6 



a Assuming that the intensity of the stellar UV radiation field in 
the DLAs is similar to the mean value measured in the 
interstellar medium of the Milky Way near the Sun. Upper-limits 
correspond to one-sigma confidence bounds. 
b One-sigma upper-limits on Wq 3935 were used in the models. 
c We use the log ./V(Crll) value given in the literature for 
modelling the DLA towards 0256+0110 as our measurement 
is highly uncertain. 

Figure [SJ). For comparison, in the lone strong Call system 
in our sample, Wq 3935 is greater only by a factor of about 
two, despite having column densities of Znll and Cm one 
order of magnitude higher. It will be of interest to verify 
the extent to which Wq 3935 and iV(H2) are related as more 
measurements of these quantities become available. 



5 SUMMARY AND CONCLUSIONS 

We have used the WHT telescope on La Palma to record at 
intermediate resolution the optical spectra of QSOs known 
to lie behind DLAs and subDLAs in the redshift range 
0.6 < Zabs < 1.3. Our principal aims were twofold: (i) to 
measure the strength of the Call AA3935, 3970 absorption 
lines in known DLAs, with a view to establishing to what 
extent DLAs and strong Can absorbers overlap; and (ii) to 
measure the Znn AA2026, 2063 and Cr n AA2056, 2062, 2066 
multiplets in the (sub)DLAs, thereby increasing the number 
of metallicity and depletion determinations in (sub) DLAs at 
low-redshift. By design, these complementary measurements 
allow us to test whether the strength of Call absorption is 
related to the metallicity of the g measured by the 

abundance of Zn, and/or to the degree of dust depletion of 
refractory elements (both Ca and Cr are readily incorpo- 
rated onto dust, whereas Zn is not). From consideration of 
such measurements in a sample of 16 DLAs and six subD- 
LAs, augmented by similar data for three DLAs and four 
subDLAs from published studies, we reach the following con- 
clusions. 

1. The new measurements presented here strengthen pre- 
vious conclusions on the metallicities and dust depletions 
in DLAs at z < 1.3. We point out that, after accounting 
for the fact that at these redshifts DLAs have mostly been 
selected from Mgn-strong systems - and may thus be bi- 
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ased in favour of metal-rich absorbers - the already weak 
redshift evolution of the mean DLA metallicity is reduced 
further: [{Zn/H}] is only a factor of about two greater at 
0.6 < z < 1.3 than at higher redshifts. We also confirm the 
dependence of dust depletion on metallicity uncovered by 
earlier studies. 

2. Most DLAs exhibit detectable CallA3935 absorption 
in our spectra. Our new data indicate that selecting absorp- 
tion systems above the threshold Wq 3935 = 0.15 A isolates 
systems with iV(Hl) > 2 x 10 20 cm -2 (although with some 
contamination from subDLAs, as evidenced by detections 
reported in the literature). However, while it seems likely 
that strong Can absorption is associated with DLAs, the 
converse is not true: approximately one third of the DLAs 
in our sample have Wq 3935 < 0.15 A, while adopting the cri- 
terion W^ 3935 > 0.30 A would miss ~ 3/4 of the DLAs. 

3. Only four out of the 19 DLAs in the combined sample 
meet the selection criterion Wq 3935 > 0.35 A of the SDSS 
survey by WHP6, and only one out of 19 would be classed as 
a 'strong' Ca II absorber by their definition (Wq 3935 > 0.5 A). 
Evidently, there is little overlap as yet between confirmed 
DLAs and the SDSS candidates identified via Can. This is 
not surprising, however, since - as demonstrated by WHP06 
- dust associated with many of the strong Call absorbers 
dims the light from the background QSOs sufficiently for 
them to be excluded from the magnitude-limited surveys 
conducted with HST up to now. 

4. Overall, we find no convincing correlation of the 
strength of Ca II absorption with the column density of neu- 
tral gas, the equivalent width of Mg II A2796, the metallic- 
ity of the gas (as measured from the abundance of Zn), or 
the degree of depletion of refractory elements (indicated by 
the [Cr/Zn] ratio). Despite the fact that, in the interstellar 
medium of the Milky Way, both Ca and Cr exhibit large 
and variable depletion factors from the gas-phase, DLAs 
with near-solar values of [Cr/Zn] do not show significantly 
stronger Call lines than those where most of the Cr is pre- 
sumably in solid form. On the other hand, the finding that 
the only 'strong' Call absorber in our sample is the most 
metal-rich DLA yet discovered suggests that there may be 
some connection between the metallicity and Call equiva- 
lent width. 

5. Presumably, a more complex mix of factors than sim- 
ply Hi column density, metallicity, and dust depletion, de- 
termines the value of Wq 3935 in DLAs. We have presented 
preliminary evidence that the 'joker in the pack' may be the 
variable fraction of Ca which is singly ionised in Hi regions, 
where Cam is often the dominant ion stage, depending on 
the densities of particles, far-UV photons, and on the tem- 
perature of the gas. 

With the data at hand, we exploit the ionisation balance 
of Ca to place constraints on the ratio of UV photons and 
particles. If we further assume that the interstellar radiation 
field in the DLAs is of the same order as that at our location 
within the Milky Way - as indeed found for DLAs at higher 
redshifts - we can arrive at estimates of, or limits on, the gas 
density n(HTOT). Our modelling, using the photoionisation 
code CLOUDY, indicates values n(HTOT) ~ 1 cm -3 for most 
DLAs with Wq 3935 £ 0.35 A. Furthermore, for an 'average' 
Can-absorber DLA (as described by WHP06) as well as 
the lone 'strong' (as defined by those authors) Ca II system 
in our sample, we deduce w(Htot) ~ 10 cm -3 . Such high 



values are similar to those which seem to be typical of H2- 
bearing DLAs at z > 2, raising the possibility that in the 
strong Can absorbers selected by WHP06 at z < 1.3 we 
may be seeing the lower redshift counterparts of the subset 
of DLAs with a measurable molecular fraction. Indeed, our 
models predict molecular fractions similar to those reported 
for relatively metal-rich DLAs at high-z, with the strongest 
Ca 11 absorbers exhibiting values consistent with the highest 
reported measurements of fu 2 in DLAs. 

If borne out by additional studies, the possible con- 
nection of strong Call absorbers to metal- and relatively 
molecule-rich gas has intriguing consequences when consid- 
ered in light of the imaging results of Hewett & Wild (2007) . 
It would mean that strong Call absorbers select gas with 
these properties not only in galactic disks, but also at dis- 
tances of tens of kpc from relatively bright galaxies. Such 
a conclusion would be highly relevant to models of galactic 
outflows and star-formation feedback. 

The currently observed population of DLAs are not, 
in general, strong Call absorbers. However, if most or all 
strong Call absorbers are DLAs and the true incidence of 
strong Can absorbers is as high as calculated by WHP06, 
these intriguing systems may represent a (potentially metal- 
rich) subset of DLA absorbers that are being missed from 
the current surveys. It will be possible to make progress 
towards clarifying this and other questions raised by the 
present work with COS measurements of metal lines and 
AT (Hi) in a sample of known strong Call absorbers. Like 
the rest of the UV-spectroscopy community, we look forward 
with anticipation to the forthcoming HST servicing mission 
that will install this much needed instrument. 
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APPENDIX A: QSO COORDINATES, 
MAGNITUDES, AND REDDENING 
ESTIMATES 

Here we list the celestial coordinates, magnitudes, and esti- 
mates of E(B — V) for the quasars listed in Tables 1 and 2. 
Magnitudes are given for the SDSS i' filter, when possible. 
E(B — V) for systems in our sample with an available SDSS 
spectrum is estimated exactly as in Wild & Hewett (2005). 
Note that the scatter in individual reddening values (due 
to the intrinsic variation in the quasar SEDs from object to 
object) cause some 'negative reddenings' to occur. 



Name 


Ha (J 2000) 


Dec (J 2000) 


m i(AB) 




- V) 


0139-0023 


01 39 38.71 


-00 23 48.1 


18.57 


-0 


099 


0253+0107 


02 53 16.46 


+01 07 59.8 


18.56 


+0 


001 


0256+0110 


02 56 07.24 


+01 10 38.7 


18.42 


-0 


031 


0449-168 


04 52 14.23 


-16 40 16.8 


m v = 17.0 






0454+039 


04 56 47.2 


+04 00 53 


mv = 16.5 






1007+0042 


10 07 15.54 


+00 42 58.4 


18.75 


+0 


029 


1009-0026 


10 09 30.47 


-00 26 19.2 


17.28 


-0 


118 


1010+0003 


10 10 18.20 


+00 03 51.2 


18.03 


-0 


039 


1107+0048 


11 07 29.04 


+00 48 11.2 


17.19 


-0 


036 


1107+0003 


11 07 36.68 


+00 03 29.6 


18.04 


+0 


089 


1220-0040 


12 20 37.01 


-00 40 32.4 


18.33 


-0 


022 


1224+0037 


12 24 14.30 


+00 37 09.0 


18.53 


-0 


052 


1225+0035 


12 25 56.62 


+00 35 35.0 


18.52 


+0 


072 


1323-0021 


13 23 23.79 


-00 21 55.3 


17.64 


+0 


220 


1727+5302 


17 27 39.03 


+53 02 29.2 


18.04 


-0 


028 


1733+5533 


17 33 22.99 


+55 33 00.9 


17.86 


-0 


018 


2149+212 


21 51 45.9 


+21 30 14 


m v = 19.0 






2353-0028 


23 53 21.63 


-00 28 40.7 


18.23 


+0 


025 


0515-4414 


05 17 07.61 


-44 10 56.2 


772b = 15.0 






0738+313 


07 41 10.70 


+31 12 00.2 


16.66 






1436-0051 


14 36 45.06 


-00 51 50.5 


18.25 






2328+0022 


23 28 20.37 


+00 22 38.2 


17.78 






2331+0038 


23 31 21.82 


+00 38 07.4 


17.51 






2335+1501 


23 35 44.19 


+15 01 18.4 


18.22 







